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Abstract 

We discuss the violation of the equality of the gauge coupling and gaugino 
coupling in the presence of soft supersymmetry breaking. Although this is 
a hard supersymmetry breaking effect, there are finite contributions to this 
difference which can be determined as a function of the soft supersymme- 
try breaking masses. The largest sources of this difference are the "super- 
oblique" corrections, which can be enhanced by a multiplicity factor and a 
logarithm of the soft supersymmetry breaking mass. This is to be contrasted 
to standard oblique corrections to the electroweak sector which decouple for 
large supersymmetry breaking. We show that these parameters can be rel- 
atively large, particularly in models of supersymmetry breaking motivated 
by solving the supersymmetric flavor problem. We also perform a detailed 
study of the non-oblique corrections for the example of squark decay. We find 
that they can be significant, and should be accounted for in the theoretical 
prediction. 



1 Introduction 



Supersymmetry has yet to be discovered. If /When it is, the information 
about the underlying theory will be limited but powerful. For example, the 
spectrum of superpartners might tell us about the mechanism of communi- 
cation of supersymmetry breaking. It is conceivable that there is very large 
supersymmetry breaking within multiplets and that there are many as-yet 
undetected multiplets carrying standard model gauge charge. Although the 
heavy superparticles will elude direct detection, there can be large virtual 
effects which could be detectable through the difference in the gauge and 
gaugino coupling. 

In this paper, we will discuss a new category of precision measurements 
which one might use to probe supersymmetric theories. As we will discuss 
in detail in Sec. |2|, these measurements have much in common with the well- 
studied "oblique" corrections which characterize corrections to the standard 
model. In supersymmetry, the gauge coupling g is equal to the gaugino 
Yukawa coupling g when supersymmetry is exact. However, when supersym- 
metry is broken, components of the supermultiplet (for example, q and q) 
no longer have the same masses, which leads to different contributions to 
gauge and gaugino propagators. This is the supersymmetric analog of the 
"oblique" corrections and gives rise to finite differences between g and g. 
The experiments will certainly be more difficult because they require very 
precise measurements and good measurements of other supersymmetric pa- 
rameters. However, these measurements could yield invaluable insights into 
the underlying physics of supersymmetry breaking. 

The key to understanding the quantities we wish to consider is to distin- 
guish hard and soft supersymmetry breaking parameters. Supersymmetry 
can only be broken softly if it is to serve as a solution to the hierarchy prob- 
lem. The focus of most research on supersymmetric phenomenology is on the 
soft supersymmetry breaking parameters, and in particular the superpartner 
masses. Hard supersymmetry breaking, on the other hand, is "forbidden". 
However, this is not true, even in the softly broken supersymmetric theory! 
The only constraint is that the coefficients of hard supersymmetry breaking 
operators are finite. There are no available counterterms in the supersymme- 
tric Lagrangian. Nonetheless, finite corrections are permissible and present. 
This is precisely analogous to the precision electroweak parameters S, T, and 
U which are parameters of the effective theory which can be calculated as a 
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function of the underlying heavy masses. In fact, our focus here is on new 
parameters which reflect hard supersymmetry breaking in the form of oblique 
corrections in the gauge sector as well. As with the precision electroweak pa- 
rameters, these parameters will serve as a test of the supersymmetric model, 
and give insights into non-standard physics. While the standard oblique cor- 
rections were sensitive to large violations of custodial SU(2) symmetry, the 
super-oblique corrections are sensitive to large splittings within a supersy- 
mmetric multiplet. 

The largest effects will occur if some superpartners have multi-TeV masses, 
though there are potentially measurable effects even with smaller mass split- 
tings. Examples of models with large mass splittings are given in |], ^ ^ [| 
in which the partners of the first two generations are heavy, suppressing po- 
tential flavor changing neutral current (FCNC) problems. The squark masses 
could also be heavier than the slepton masses because of the larger contri- 
bution from the running of the masses due to the SU(3) strong coupling 
constant, or because of a larger tree- level mass in a gauge- mediated model. 
In Sec. |]we present three different models which give different predictions of 
the super-oblique corrections motivated by these considerations and compare 
the effects. 

Oblique corrections are also important because they are process indepen- 
dent, and therefore can be constrained by different experiments. Non-oblique 
corrections on the other hand depend on the particular process under con- 
sideration. Furthermore, virtual heavy states decouple and will not give 
large effects. Finally, there is no multiplicity factor (if there are many mul- 
tiplets with large splittings). However, it is possible that non-oblique super- 
symmetric loop corrections could be important when extracting the oblique 
parameters. This is because of the enhancement by a group theory factor 
(from diagrams with internal gauge bosons and gauginos). Since the domi- 
nance of the oblique corrections requires that these effects are small, we have 
also calculated non-oblique corrections. These include vertex, wave function 
renormalizations, and gluon emission. A detailed check on our calculation 
is the cancellation of infrared divergences, as is discussed in Appendix [EJ 
For simplicity, in this paper, we apply our calculation of non-oblique correc- 
tions only to squark decay; we find the non-oblique corrections can in fact 
be significant. This indicates that it will be important to also incorporate 
non-oblique corrections in order to extract the masses of heavy parameters 
dominating the oblique-corrections. 
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Our work is in some ways complementary to the recent papers by Cheng, 
Feng, and Polonsky|| and Pierce, Nojiri, and Yamada|| who considered in 
greater detail methods for extracting oblique parameters from experiments. 
We instead focus on the theoretical "uncertainty" arising from non-oblique 
corrections which could corrupt the extraction of parameters if they are large. 
Of course, these can ultimately all be incorporated. We also add to the con- 
sideration of experimental signatures by giving the prediction for the produc- 
tion of bbg and ttg (assuming it is kinematically accessible) which could prove 
to be a promising method of extracting SU(3) gaugino Yukawa coupling if 
the supersymmetric partners are light. 

We begin in Sec. || discussing the analogy between the super-oblique and 
oblique corrections and the non-decoupling effect of heavy superparticles. In 
Sec. ||, we present different examples of models that contain heavy super- 
partners. We give the results for super-oblique correction for these models 
in Sec. f|. In Sec. [5], we consider additional gauge representations and their 
contribution to the super-oblique corrections. In Sec. |, we discuss possible 
ways to measure the super-oblique parameters and give the results for the bbg 
and ttg production cross sections. In Sec. [7], we investigate the non-oblique 
corrections and show that they could be quite significant compared to the 
oblique corrections. In Sec. |8|, we conclude. We have several appendices 
with details of our calculations. Appendix A contains the 0(a) corrections 
to vertices and self-energies for arbitrary parameters. Appendix B gives the 
super-oblique corrections. Appendix C gives the results of gluon emission, 
which will be necessary to cancel infrared divergences. In Appendix D, we 
consider the results in the limit of exact supersymmetry. Appendix E shows 
that the results of Appendices A and C combine to give infrared finite re- 
sults. The final appendix gives the results for super-oblique corrections from 
a messenger sector. 

2 Oblique and Super-Oblique 

It is useful first to review standard oblique corrections. Precision electroweak 
corrections were important in that they tested consistency of the standard 
model and gave insight into the high-energy world. The oblique corrections 
are the corrections to the gauge boson propagator. The parameters which 
characterize the most important effects can be obtained by retaining the 
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leading pieces in a derivative expansion, accounting for gauge invariance. The 
six "parameters" are II' , II' z , llVjy, njy^, Hzz, tl'zz- We chose to work in 
the basis of tree-level mass eigenstates and retained only the nondecoupling 
pieces of the propagator, those of dimension four or less. It is useful to 
absorb three of the parameters into g, g', v, which leaves three new finite 
parameters which are forbidden by gauge invariance but are permitted once 
it is spontaneously broken. In Ref . , these three parameters were identified 
as as SZzz, fiZ^z, Smzz- In the more widely utilized Peskin-Takeuchi naming 
convention, they were absorbed into parameters [§, || S, T, U, while a third 
popular convention from a paper of Altarelli and Barbieri [|10[ absorbed them 
in the three parameters ex, €2, 63. The importance of these parameters is 
that even heavy states do not decouple. For example, one can compute the 
contributions of a heavy doublet to the standard S parameter, 



Of course, the many new particles present in a supersymmetric theory can 
contribute to the standard oblique corrections. In general, these effects are 
small because heavy superpartners decouple. However, there is a potentially 
more interesting class of "oblique" corrections for supersymmetric theories. 
Supersymmetry guarantees equality of the gauge and gaugino Yukawa cou- 
plings. When supersymmetry is broken, these are no longer guaranteed to 
be equal at low energy. 

We consider the relative difference of the gauge and gaugino Yukawa cou- 
plings, which is sensitive to supersymmetry breaking. Due to supersymmetry, 
this difference is UV-finite and hence independent of the renormalization scale 
fi. This difference could be measured by comparing the Yukawa coupling g, 
extracted from some physical process at energy scale q 2 , to the gauge cou- 
pling g at that same energy. Because a is reasonably small, and logarithms 
are not large, we choose to work to leading order in a, rather than taking an 
effective theory approach and performing a full leading log calculation. This 
is useful since the finite pieces can be competitive with the logarithmically en- 
hanced terms, while higher order terms are negligible. Note also that these 
parameters are //-independent and well defined, so one can calculate both 
logarithmically enhanced and nonlogarithmically enhanced contributions. 

We define the gauge coupling g(q 2 ), and the gaugino Yukawa coupling 
g{q 2 ), as the couplings containing all the "super-oblique" corrections; in 




(2.1) 
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other words, as the couplings containing contributions to the gauge (or gaug- 
ino) propagator from all possible squarks and quarks running in the loops. 
We note that our definition of g(q 2 ) (which is process independent) differs 
from the convention of Ref. ||11|| , who used an effective theory language and 
therefore did not include in their definition oblique contributions from light 
fermion loops. Our two conventions differ by these finite contributions, which 
are generally smaller than the logarithmic contribution but not by a substan- 
tial amount. Since one can readily calculate these additional loops, which 
the full order a calculation would require in any case, we choose to include 
these pieces also in our definition of oblique corrections. 

Hence, in our case the difference between the couplings will receive contri- 
butions from all "super-oblique" diagrams (including the light fermion loops). 
These consist of the quark and squark contributions to the the gauge prop- 
agator (Fig. |6| (3) and (4)), and the gaugino propagator (Fig. [5] (2)). For 
general choices of superpartner masses, the difference at one loop takes the 
form of a function 

g(q 2 )-g(q 2 ) « (K K ™j\ a F g3 + F gi -F- g2 

g(q 2 ) 4vr ° b ^ q 2 ' q 2 ' q 2 ) 4vr 2 ' 1 ] 

which depends on the ratios of the various masses with the scale q 2 (F g s, 
FgAi Fg2 are explicitly given in Appendix |A| as integrals over Feynman pa- 
rameters). For example, in the case of the SU(3) gauge and gaugino Yukawa 
couplings, the contribution of a single heavy doublet superfield to this differ- 
ence at q 2 ~ m~ is 

93 - 93 «3 / T _ _ / m g 



93 24vr 



Log — -T7T • (2-3) 



This formula is only true in the limit m q 2 <C q 2 and q 2 <C nig 2 . A detailed 
derivation of this formula is given in Appendix [B[ Also notice that because 
physical quarks are generally lighter than the gauginos, the logarithm is cut 
off at the gaugino mass. Therefore, the logarithm is only large if there is a 
sizable splitting between the squark or slepton and gaugino masses. 

Although the oblique corrections first arise at one loop (and there is an 
additional 1/6 suppression compared to the usual a/(47r)), there are factors 
which can make it large. First of all, we see that there is a logarithm; if 
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there is a large mass splitting, there will be a corresponding enhancement. 
Second, there is a multiplicity factor; "every" particle with gauge charge 
contributes. The logarithm here differs from standard oblique corrections 
in that decoupling works differently. Heavy scalar superpartners contribute 
more rather than less. 

If we call the soft scalar mass ttiq and the coefficient of the trilinear scalar 
coupling A, for the standard oblique corrections (e.g. the S parameter), it is 
easy to see that 

AS oc — 2 > ( 2 - 4 ) 

m z rriQ 

which decouples for large m . 

In the case of the super-oblique corrections which we discuss, one finds 

2 

6g(m q ) oc Log— §, (2.5) 

(where we now assume the quark mass is bigger than the gaugino's and take 
q 2 ~ m q 2 ) which clearly does not decouple, and is in fact logarithmically 
enhanced. This logarithm can however be readily understood ||12|| . If we 
assume a chiral multiplet in which the scalar partner is heavy, there is an 
energy regime between the fermion (or gaugino) mass and the heavy scalar 
mass in which the gauge and gaugino coupling run differently. In fact one can 
also determine that when the scalar is the heavy partner, the gauge coupling 
at low energy will be reduced relative to the gaugino couplings due to the 
sign of the fermion contribution to the (3 function for the gauge coupling. 
Decoupling does occur however for super-oblique corrections when m q ~ 



rrig 



(m 2 q — mjp 



5g ex v q - q ' . (2.6) 



m q 

Notice here that we have assumed the squark and quark to be nearly degener- 
ate and heavier than the gaugino. Therefore what is relevant to super-oblique 
corrections are standard model gauge charge multiplets with large supersym- 
metry breaking splitting in the masses. If these mass splittings are not very 
large, the constant term 11/12 in Eq. |2.3| can be somewhat significant. For 
instance, if m q = 5m,g the constant term is a 30% correction to the logarithm. 

We conclude that heavy scalar partners can lead to interesting deviations 
from the exact supersymmetric prediction g = g. 
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3 Large Effects and FCNC 



It is reasonable to ask whether there is any reason to believe that there 
are large mass splittings. The answer might be yes, given the current ideas 
for resolving the FCNC problems in supersymmetric theories. The apparent 
degeneracy of squarks or alignment with quarks is still not understood. There 
are two types of scenarios which offer the possibility of large mass splittings 
among the superpartners. It could be that squarks (and gluinos) are heavier 
because of the larger SU(3) coupling. Alternatively, a suggested resolution 
of the FCNC problem involves heavy first two generations [13], 0, §. 

An example of a class of models in which the colored states are heavy is 
the gauge-mediated models which have been suggested as a resolution of the 
FCNC problem. It could also be that colored scalars get enhanced masses 
from renormalization group running. Of course, either of these scenarios 
would probably make both the squarks and gluino heavy. A very heavy 
gluino would not be subject to the experimental tests we are discussing, 
although there could still be measurable oblique corrections for the SU(2) 
and U(l) sectors. Just to test how big the effects can reasonably be expected 
to be in a model in which the differences between the gauge coupling and 
gaugino Yukawa coupling can all be measured, we took a toy model with 
heavy squarks and a relatively light gluino. Although not motivated by 
any particular model, it serves to illustrate the potential size of the oblique 
corrections. 

There are several groups who investigated the possibility of heavy first 
two generations. An example of such a model is [0, where an anomalous 
U(l) symmetry which served as a mediator of supersymmetry breaking was 
introduced. If only the first two generations are charged, their mass can 
be of order Air bigger than the lighter third generation states. The gaugi- 
nos are also among the light states. In Ref. ||, the heaviness of the first 
two generations is achieved through the constraints of naturalness with some 
restrictions on fine tuning. All the superparticles can be divided into two 
groups. The first, called "the brothers of the Higgs", contains superpartners 
whose masses will strongly influence the mass running of the Higgs and thus 
are strong constrained by naturalness. The second, called "the cousins of the 
Higgs", consists of superparticels which are less constrained and hence can 
be heavy. The first two generation squarks are in the latter group and their 
mass can be in the range of 900 GeV— 5 TeV. Related ideas are incorporated 
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into the "More Minimal Supersymmetric Standard Model" (also called effec- 
tive supersymmetry) |IJ , in which one attempts to allow the maximal masses 
consistent with naturalness bounds. In this model, a new gauge group G is 
introduced, which enlarges the accidental symmetry group and thus forbids 
the renormalizable B and L violating interactions. It also introduces a new 
mass scale M ~ 5 — 20 TeV which sets the mass scale of the first two genera- 
tions. This model permits the first two generations of squarks to be heavier 
which goes a considerable way towards solving flavor changing problems. In 
addition, the requirement of naturalness implies that some squarks (tj,, &l) 
and most gauginos must have a mass below ~ 1 TeV. One can then imagine 
mi )2 ~ 20 TeV, m 3 ~ 100 GeV - 1 TeV. 



4 Parameters and Numbers 

One can now ask for the magnitude of the finite difference between the gauge 
and Yukawa couplings for the kinds of models described above. There are 
three measurable quantities associated with g 3 , g 2 , gi, which are the coupling 
constants of the gauge groups SU(3), SU(2) and U(l) respectively. These fit 
in well with a modified Peskin-Takeuchi naming convention v = (gs — gzj/gz, 

w = (g 2 - g2)/9* y = (gi - 9i)/gi- 

We now present the magnitude of each of these quantities for the types 
of models described in Sec. || in which there are mass splittings. We make 
simple assumptions below to get an idea of the magnitude of the effects. In 
the following we will use the value of the various couplings at Mz to estimate 
the effects, and we take q 2 ~ m g a ugmo- 

We define three different models: 

Model 1 Only the squarks are heavy. We take = f^wl m q = 0-05, and 
since this is a gauge mediated model nig : = a 3 : a 2 « 7 : 2, where 
a 3 = 0.12, a e = 1/128, sin 2 W = 0.23. Thus, R\ = m~gjm~ q = 0.18. In 
Table [I], we present the numerical results for Model 1. Here we have 
ignored the additional finite contribution from the sleptons. If, for ex- 
ample, mj = rriyy, one calculates the finite contribution a2/247r(ll/2 — 
VSti/2) instead of a 2 /247r(-ll/12) in Eq. |J (where we consider SU(2) 
case, and similarly for U(l) case). 

Model 2 A More Minimal Spectrum, in which the masses of 6r, u, d, s, 
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~v 3.82% log R g 9 - 1.8% -8.4% 
w 0.81% log Rf -0.37% -2.8% 
_y 0.30% log Rf -0.14% -1.0% 



Table 1: Super-oblique corrections for Model 1. 



v 0.95%log J R[ + 2.86%logi?|-1.8% -14.3% 

w 0.36%lo gj Rf + 0.72% log Rj - 0.5% -5.1% 
y 0.16%logfl{f + 0.38% log gf -0.25% -2.6% 

Table 2: Super-oblique corrections for Model 2. 



c, e, /i are mn = 20 TeV, whereas the masses of &l, r are tnt =1 
TeV, while = lOOGeV, and m g = 350GeV. In this model, R[ = 
iTig/mi, = 0.35, R^i = rrig/mn = 0.018, = m vy/ mL = ^-1 and 
R¥( = m^/ffln = 0.005. In Table |2], we present the numerical results 
for Model 2. 

Model 3 A modified more minimal model similar to Model 2, in which the 
sleptons have mass =1 TeV. In Table we present the numerical 
results for Model 3. 

We can see that v is the same for Models 2 and 3 since changing the slepton 
mass would not affect the SU(3) coupling constant, which only gets contribu- 
tions from q, q and g. v is much larger than w and y since the corresponding 
coupling is much larger. We write the super-oblique corrections in terms of 
the ratio of the mass splittings in Tables [l], 0, |3], where it is evident that 
the finite piece 11/12 does give a sizable contribution if the logarithm is not 
much larger than 1. The results we get for Model 1 are consistent with the 
results for the "heavy QCD models" in Ref. [0 when the gluino is light. 



Our model 2 is close to the "2-1 models" [14| except that we put 6r in the 



0.95% log Rl + 2.86% log R\ - 1.8% -14.3 % 
0.54% log Rf + 0.54% log Rf - 0.5% -4.6% 
0.32% log Rf_ + 0.22% log fljf - 0.25% -2.2% 

Table 3: Super-oblique corrections for Model 3. 
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heavy sector. When we let 6r also be light, the coefficient before logRjj in 
Model 2 also fits well with the results given in Ref. [14 . However, 



we also calculate the contribution from the finite piece 11/12 and the light 
squarks, whose masses are still large enough (compared to the gauginos) to 
give noticeable contribution. 



5 More Non- Standard Physics 

In this section, we consider the role of the super-oblique correction measure- 
ments not only as a test of the parameters of the minimal extension of the 
standard model but also as a potential probe of other new physics. We have 
already stressed the fact that heavy squarks can be large. This is especially 
important as it is precisely when the squarks become kinematically inaccessi- 
ble that their virtual effects become largest. If one does not find the squarks, 
or if they are at the edge of the kinematically accessible regime, one would 
like to be able to confirm their existence and attempt some estimate of their 
mass through their virtual effects. 

However, it might be that the squarks are relatively light, and one still 
measures large deviations from the supersymmetric relations. Provided they 
are not too large, one might ask how this could be consistent with super- 
symmetry. In fact, this is probably the most exciting possibility. It would 
indicate the existence of nonstandard fields which participate in supersym- 
metry breaking and therefore have a large mass splitting relative to the mass 
of the fermions or scalars and which carry standard model gauge charge. 
An obvious candidate in terms of existing models is a messenger sector, 
although possibilities extend beyond this. As an example, suppose the mes- 
senger sector consisted of a large vectorlike representation under SU(5). One 
should note that for the "standard" messenger scenario, with a singlet cou- 
pled to messengers through a Yukawa term, the masses of the scalars are 
\?\ 2 S 2 ± \F S whereas the fermion mass is XS. In this case, q 2 (the exter- 
nal four-momentum squared for g and g) is taken to be smaller than the 
fermion mass. We take the limit q 2 <C ^f erm i on 5 <? 2 ^ m scaiar to calculate 
the super-oblique corrections. Let m 2 cal ^ T /m 2 crmion — 1 ± x. The effects are 
unfortunately small unless x is very close to 1 (see Appendix |F[). Without 
an enormous multiplicity factor, a nontrivial mass splitting is needed to get 
a sizable contribution. 
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6 Measurements of Oblique Parameters 



It is clear that there can be fairly sizable effects. It is therefore extremely in- 
teresting to ask the question whether the parameters can be measured. The 
numbers presented above serve as benchmarks for interesting measurements, 
although a measurement at any level of accuracy will constrain our confi- 
dence in supersymmetry or equivalently constrain extensions of the minimal 
scenario. 

Several groups have begun to investigate this question. Initially this was 
pursued as a test of supersymmetry at tree level Jll| and it was later studied 
in order to extract squark masses if they are heavy |T2| , [J. More recently, 
several groups, |L5|, [U], 0, @] have studied super-oblique corrections. We 
summarize some suggestions for measurements. We also present new results 
for the bbg and ttg productions. 

One measurement which has been suggested of the U(l) coupling is in a 
paper by Nojiri, Fujii, and Tsukamoto where they studied cr(e + e _ — > 
epep). This process proceeds through s-channel gauge exchange and t- 
channel neutralino exchange. One can measure gi and M\ by measuring the 
angular distribution of the differential cross section da(e + e~ — > e^e^) / d cos 9. 
It has also been suggested by Cheng, Feng, and Polonsky M to look at the 
e~e~ scattering cross section through t-channel neutralino exchange when 
running the collider in the appropriate mode. This can lead to even greater 
precision since some of the troublesome backgrounds are absent and the 
ability to highly polarize the beams increases the srSr cross section further 
[ O , It should be borne in mind that although the U(l) coupling might 
be the most accurately measured, the magnitude of y is the smallest of the 
three precision measurements. It is therefore worthwhile to investigate the 
possibility of measuring w and v as well. 

Several ways of measuring w have been suggested; the most appropri- 
ate will depend on the region of supersymmetry parameter space which is 
realized O, |5| . We summarize some of the suggestions of Ref. 0] . If the 
charginos are kinematically accessible, their production cross section can be 
used to measure w, for example if they are pair produced with a contri- 
bution from sneutrino exchange. Neutralino production (and decay if more 
than one mode is sizable) can also be used to test w (and y). It might also 
be possible to test the gaugino coupling through the chargino mass depen- 
dence on the HHW vertex, though this will help only away from the region 
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where charginos are nearly pure gaugino. It should be noted that a good 
measurement of will be important when constraining w to high preci- 
sion by some of these methods. Another possible way to study w is through 
the branching fraction of charginos if there are appropriate modes which are 
kinematically accessible. Finally, selectron and sneutrino production could 
provide important probes. 

Finally we consider the measurement of v , which requires measuring the 
SU(3) Yukawa coupling. It is very likely that the best test will be at a hadron 
collider, where gluino and squark can be produced abundantly. In the case 
that squark decay is observed, the branching ratios into gluinos and winos 
can serve as probes of the super-oblique corrections. Notice that the gluino 
production section goes primarily through the gluon coupling, and therefore 
does not serve as a test at all. The branching fraction tests are also possible 
at an e + e~ collider if squarks are sufficiently light. 

At an e + e~ collider, it might be possible to test the 2-1 type models if t, g 
are sufficiently light@ or if b, g are sufficiently lighti. We consider the the cross 
section for ttg and bbg production to get an estimate of the possible sensitivity 
to v . For b quark and squark production one predicts a few hundred events 
according to the gluino and squark mass (below a few hundred GeV) and 
roughly the same for the top quark. In Fig. |8] and Fig. |9] we present the 
numbers of events produced per year at an e + e~ collider at E cm = 1000 GeV 
and luminosity L = 5 x 10 33 cm~ 2 s -1 . The number of events decreases with 
the increase of mq and rrig due to the decrease in phase space. The number of 
b events is slightly less than t events, though their phase space is larger. This 
is because of the different couplings of the b and t to Z and 7. When mq > rrig, 
it is possible that the intermediate squark state is on shell, which gives a peak 
in the production rate which we have not shown. This is more significant to 
bbg production since ~ 4 GeV, so there is a larger range of parameters 
for which the intermediate b is on shell. In this case one is in fact measuring 
the squark decay branching ratio into gluino, which can also determine the 
gaugino Yukawa coupling. By measuring the branching ratio for the decay 
of the q into q — > qg, q — > q'x* 1 and q — > qx°, one can hope to extract the 
coupling, if the other parameters (masses) can be independently determined 
0. The number of events we present here is the event production, not 

a We thank Michael Peskin for suggesting this possibility 

b We thank Jonathan Feng for pointing out this alternative possibility. 
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corrected for cuts necessary to remove background events. Nonetheless for 
sufficiently light gaugino and third generation squarks, it seems likely that 
a measurement of the gaugino coupling at at least the 15% level should be 
possible. 

7 Non-oblique Contributions 

When extracting the qqg coupling from a physical process, such as squark 
decay branching fractions, it is important to consider finite contributions 
from non-oblique diagrams. Although one usually ignores such terms, they 
could be significant because of the C 2 (G) factor (C 2 (G) = N for SU(N)). 
For instance, in the non-abelian gluon and gluino loop diagrams (Fig. ^| (4)), 
even though there are no large logarithms, the C^G) acts as a multiplicity 
factor, enhancing their effect. Hence, it is useful to check the relative size of 
the non-oblique corrections. 

It is important to bear in mind that the /z dependence in the oblique pa- 
rameters cancel; everything is finite. We will sometimes refer to finite pieces; 
by this we mean the nonlogarithmically enhanced contribution (that is, the 
finite pieces in an effective theory approach). The finiteness is important in 
that it means that the answer at one-loop is renormalization scheme inde- 
pendent. The dependence on momentum is through the dependence on the 
physical momenta which enter the loops. For any given process these are 
specified. For example, in the case of squark decay, the mass squared of the 
squark will enter. In a process with a virtual squark, the loops will be func- 
tions of q 2 and the squark mass squared; the q 2 dependence of the virtual 
corrections can be incorporated when one integrates over q 2 . Because it is 
clearly simpler, we illustrate the non-oblique corrections (which are process 
dependent) with the specific case of squark decay. One can however apply 
our results to any process by which the oblique parameters will be measured. 

Now let us focus on the q decay process q — > qg. The relevant non- 
oblique diagrams include the quark and squark self energies (Fig. [|, Fig. |5|), 
the non-abelian gluon and gluino self energies (Fig. || (1), (2) and Fig. |3] 
(1)), the qqg vertex diagrams (Fig. ^|), and the gluon emission diagrams 
(Fig. 0). Here, all the diagrams are necessary to render a gauge invariant 
answer free of infrared divergencies. The quark and gluino are assumed to be 
always on shell, while the squark could be offshell at energy q 2 . We note that 
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for completeness we have included both abelian diagrams (Fig. |^ (1), (2), 
Fig. |], Fig. H) (proportional to C^iV)), and non-abelian diagrams (Fig. |2] 
(3), (4), Fig. § (1), Fig. | (1), (2)) (proportional to C 2 (G)). We use the 
dimensional reduction regularization scheme (DR), rather than dimensional 
regularization, to avoid obtaining a finite difference between the gauge and 
the Yukawa couplings in the supersymmetric limit. This finite difference 
occurs in dimensional regularization because of the mismatch of fermionic 
and bosonic degrees of freedom for dimension d > 4. 

In addition, we introduce a gluon mass, A, to regulate infrared divergen- 
cies. Since the quark mass is usualy much smaller than that of the squark 
and gluino, we could have set m q to zero. However, we found that because 
m q effectively regulates the collinear divergences in the angular integral of 
the real gluon emission, keeping it non-zero made it easier to keep track of 
the various sources of divergences. As a result, separation of soft emission 
divergences, proportional to the logarithm of A, from those of collinearity 
(proportional to the logarithm of m q ) was straightforward. Identifying the 
infrared divergence in a given vertex diagram with the corresponding diver- 
gence in the emission diagram, we were able to explicitly verify that the sum 
of diagrams is infrared finite (see Appendix The total UV- finite virtual 
non-oblique correction is defined to be 

Fvi rt uai(g 2 , ml ml ml A 2 ) = Y,^ + \h F * + \h F v + \ F ^ ( 7 - X ) 

i=l i=l i=l 

Here, the Fy; are the vertex diagrams contributions, and F qi , F q i are 
the quark, squark, and gluino self-energies. All the Fs are explicitly given in 
terms of integrals over Feynman parameters in Appendix We also define 
the UV-finite contribution arising from gluon emission, -F emission , by 

T q _ q ~ gg (q 2 ) = £;f (g 2 )F cmission (g 2 , m 2 , m?, mj, A 2 ). (7.2) 

Here, To is the tree level decay rate for q — > qg and its explicit formula is given 
in Eq. |C.6| . Note that in the three-body final state decay we are integrating 
over all of phase space, effectively including soft and hard gluon emission. 
This constitutes an overestimate of the emission finite contribution. However, 
the sign is the opposite to the overall sign of the non-oblique corrections, so 
this constitutes an underestimate of their size. 



14 



The total squark decay width is thus given by 



Ttotal — To ( 1 + ~ (2-FVirtual + -^emission) ) • ' ' - >i 



Recall that in Eq. |2.2| 



then 



g{q 2 ) - i^ ob ' (7 ' 4) 

a = a(l-^2F oh ). (7.5) 

Combining Eq. ^3| and |775| (We use a in stead of a in Eq. ^3| because the 
difference is high order.), the squark decay width is given in terms of g(q 2 ) 
as 

/ a 

Ttotal 1^0 ( 1 ~. (2-Fvirtual -^emission 2_Fo b 

\ 47T 

= r (l + £(2F aon _ ob - 2F ob )V (7.6) 

where F non „ ob is defined by F non _ ob = F virtual + l/2F emission and T is similar 
to T except that the coupling constant is g instead of g. 

The numerical results for SU(3) (tt/47r)F non _ ob for on shell q decay is 
shown in Fig. [1^. For a range of m g /mq of 0.5 — 0.8, which is a reason- 
able range for the parameters which might be accessible to an e + e~ collider, 
(a/47r)F non _ ob changes from -8.5% to -1.57%. This can be a significant frac- 
tion of v. Furthermore, it has the same sign to -F ob contributions, thereby 
lessening the observed effect (Notice the relative minus sign between F ob and 
-Pnon-ob)- This cancellation is more significant when m g /mq is smaller. Recall 
that we are imagining a model like Model 2. For m g /m^ = 0.5 — 0.8, with 
m L /m H = 1/20, we obtain v = —10.2% — > —9.7% (where the parameter is 
given at q 2 = ml). 

It might seem surprising that the non-oblique corrections are larger in 
magnitude for smaller m g jm^ which seems to contradict the decoupling 
effect of the heavy sectors. Actually this is due to some large negative total 
shift of the non-oblique corrections with respect to zero, while .F nc , n _ ob still 
goes proportional to m g /mq, especially when the ratio is small, which shows 
the decoupling effect of the heavy sector. 
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-Pnon-ob for SU(2) can be obtained similarly, and is -1.43% for m^/mj = 
0.5 and -0.24% for m^/mj = 0.8 (as compared to w = -2.6% -> -2.4%). 
A similar calculation for only the abelian part is done in Ref. [|TB|, which is 
consistent with our result except that in their Eq. (13), the constant in F rea i 
should be 7/2 instead of 13/4 and the Fg in their Eq. (19) should have the 
opposite sign. 

We only calculated the non-oblique corrections for on-shell squark decay. 
For other physical processes, the non-oblique corrections can be done simi- 
larly using the same method. For example, in the process e + e~ — > W + W~ 
through t-channel v exchange, one would also want to include non-oblique 
corrections. However, there are additional radiation diagrams which in gen- 
eral must be incorporated. 

As an additional check of our results, we have verified that in the SUSY 
limit (m g = mg,m g = 0), the non-oblique UV-finite corrections to g\ and g- x 
are identical (See Appendix [B]). 

8 Conclusions 

To conclude, it is obvious that if supersymmetry is discovered, we would like 
to gain as many handles on the underlying supersymmetric theory as can 
be accessible. Super-oblique corrections provide a different perspective into 
the physics of supersymmetry breaking and into the consistency of the stan- 
dard supersymmetric sector. As with the standard oblique corrections, they 
test consistency of the simplest version of the theory and provide access to 
nonstandard physics. They are most useful precisely when they are largest; 
namely there are heavy superpartners which cannot be directly observed. 
The particular parameters we discussed are the supersymmetric analogs of 
the standard-oblique parameters. The one-loop gluon and gluino propaga- 
tors give contributions to the gauge and gaugino Yukawa coupling difference, 
which are finite. In nonstandard models, especially those addressing the 
supersymmetric flavor problem, it is likely that the parameters we have de- 
scribed could be large. We have shown that non-oblique corrections are also 
potentially significant, and should be incorporated in the theoretical predic- 
tions. 

The super-oblique parameters will be difficult to measure, although there 
are many potential candidate measurements. Clearly, the better measured 
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are the super-oblique parameters, the more constrained will be the under- 
lying model. These parameters merit further study, both in understanding 
experiments and in further calculations of non-oblique effects. 
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A 0(a) corrections to Vertices and Self En- 
ergies 

In Fig. p] and [| we show the one-loop Feynman diagrams that contribute to 
the corresponding gqq and qqg vertex corrections (denoted by F Vi and Fy{). 
Using the DR renormalization scheme with D = 4 — 2e and regularizing the 
infrared divergences by an infinitesimal gluon mass A, we geifi 



'VI 



(C 2 {N) - 1/2C 2 (G)) 



- + log/! 2 - 1- 



J dxdydz 5(1 — x — y — z)(2 log((x 2 + 2xy + y 2 )m g 2 — xyq 2 + z\ 2 ) + 
2(2 - 2x - 2y - x 2 - 2xy - y 2 )m 2 + 2(— 1 + x + y - xy)q 2 



x 2 + 2xy + y 2 )m q 2 — xyq 2 + z\ 2 

- + log fi 2 — J dxdydz 5(1 — x — y — z) 
2 log((x + y)rriq 2 + zrrig 2 + (—x — y + x 2 + 2xy + y 2 )m 2 — xyq 2 ) — 



(A.l) 



V2 



(C 2 (N) - 1/2C 2 (G)) 



4(x + y)(l — x — y)m q 



(x + y)rriq 2 + zrrig 2 + (— x — y + x 2 + 2xy + y 2 )m q 2 — xyq 1 



(A.2) 



V3 



C 2 (G) 



3 

2? 



log /i 2 — - — J dxdydz 5(1 — x — y — z) 



c In all the formulas and figures, we use SU(3) as an example. For SU(2) and U(l), the 
calculation is similar. 
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(3 log((— x — y + z + x 2 + 2xy + y 2 )m q 2 — xyq 2 + (x + y)A 2 ) + 
3(1 — 2x — 2y + x 2 + 2xy + y 2 )m 2 + (1 + x + y — xy)q 2 , 



(— x — y + z + x 2 + 2xy + y 2 )m q 



xyq* 



(x + y)A 2 



(A.3) 



'V4 



C 2 (G) 



1 

2i 



1 1 /• 

+ - log /i 2 — - — / dxdydz 5(1 — x — y — z) 



(\og(zm q 2 + (x + y)m g 2 - (x + y - x 2 - 2xy - y 2 )m 2 - xyq 2 ) 



nig 2 + (1 — 2x — 2y + x 2 + 2xy + y 2 )m 2 + xyg 2 
zm^ 2 + (x + y)rrig 2 — (x + y — x 2 — 2xy — y 2 )m q 2 — xyq- 



(A.4) 



Fvi = (C 2 (N) - 1/2C 2 (G)) 



1 , 2 1 

i + -2 



dxdydz 5(1 — x — y — z) 



(2 log(ymq 2 — yzm-g 2 + (z — xz)m q 2 — xyq 2 + x\ 2 ) + 

(— 2x — 2y + xy + y 2 )mg 2 + (1 + x 2 + xy)m 2 + (2x — xy)g 2 + 2zm g m Ci 



ym q 2 — yzm-g 2 + (z — xz)m q 2 — xyq 2 + xA 2 



Fv2 = (C 2 (N) - 1/2C 2 (G)) / dxdydz 5(1 - x - y - z) 

2((1 - y)mg 2 + (1 - x)m q 2 + (x + y)m q m q ) 
zm q 2 + (x — y + xy + y 2 )m g 2 + (— x + y + x 2 + xy)m 2 — xyq 2 



(A.5) 



(A.6) 



V3 



C 2 (G) 



1 

2i 



1 If 

H — log /i 2 — - — / dxdydz 5(1 — x — y — z) 



(log(xrriq 2 + (-y + z + xy + y 2 )rrig 2 + (-x + x 2 + xy)m 2 - xyq 2 + yA 2 ) + 



(1 + xy + y 2 )m g 2 + (-2x - 2y + x 2 + xy)m 2 + y(2 - x)g 2 + 2zm g m q 
2(xm q 2 + (— y + z + xy + y 2 )m g 2 + (— x + x 2 + xy)m 2 — xyq 2 + yA 2 ) 

(A.7) 

2 /" 

- + 2 log/! 2 — 1 — y dxdydz 5(1 — x — y — z) 



^V4 



C 2 (G) 



(4 log((xy + y 2 )rrig 2 + (x 2 + xy)m 2 — xyq 2 + zA 2 ) + 



(1 — x — 2yz)m g 2 + (1 — y — 2xz) 



(z + 2xy)q 2 — (x + y)m g m q 



(xy + y 2 )m g 2 + (x 2 + xy)m q 2 — xyq 2 + zA 2 



(A.8) 



Here 1/e = 1/e — 7e + log47r, /i is the arbitrary renormalization scale, x, y, 
z are the Feynman parameters, g 2 is the external four-momentum squared 
for g or q and all the other external particles (q and g) are on shell. C 2 (R) is 
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defined by t a (R)t a (R) = C 2 (R) ■ 1, where t a {R) is the generator of SU(N) in 
the representation R with normalization ti[t a (N)t b (N)} = l/2S ab = C(N)S ab 
for the fundamental representation N. 

The one-loop corrections to the g self energies are shown in Fig. 13, giving 



si 



C 2 (G) 



-- ~ log/i 2 + 



2 J dxdy 5(1 — x — y)((l — x) \og{xrrig 2 + {—x + x 2 )q 2 + y\ 2 ) 



— 4(— x + x 2 )m g 2 + 2(1 — x)(— x + x 2 )g 2 . 
xm^ 2 + (— x + x 2 )g 2 + yA 2 

F~ g2 = C(N) [-^-21og/x 2 + 

4 y dxdy 5(1 — x — y)((l — x) log(t/m 9 - 2 + xm g 
2(1 -x)(-x + x 2 )g 2 



(A.9) 



'-x + x )q ) + 



yrriq 2 + xm g 2 + (— x + x 2 )q 2 ' 



(A.10) 



Here again the g 2 is the external four-momentum squared of the g, and is the 
external four- momentum squared of the q, q or g in the q, q or g self energies 
below. 

With the same conventions we use above, we can get the corrections to 
the self energies of q, gand g: (the corresponding one-loop Feynman diagrams 
are Fig. || Fig. [| and Fig. ||) 



C 2 (N) 



1 



logyU 2 + 



2 J dxdy (5(1 — x — y)((l — x) log(xm g 2 + (— x + x 2 )g 2 + y\ 2 ) + 
2(-x + x 2 )((l-x)g 2 



^ 2 - 2m, 2 ) 



?2 



xm g 2 + (—x + x 2 )g 2 + yA 2 
C 2 (iV) r 1 



(A.ll) 



-- - logyU 2 + 

. e 

2 y dxdy (5(1 — x — y)(x log(xmqi 2 + ym g 2 + (— x + x 2 )q 2 ) + 
2x(— x + x 2 )g 2 



xrriq 1 + 2/m§ 2 + (—x + x 2 )q 2 ' 



(A.12) 
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2 „ 2 

-^-21og/i 2 + -- 
e 3 



% = C 2 (N) 

4 J dxdy 5(1 — x — y)(3(—x + x 2 ) \og(xm g 2 + {—x + x 2 )q 2 + ym g 2 ) — 
x(l — x)(— x + x 2 )q 2 



xm q 2 + (—x + x 2 )q 2 + yrrig 2 



(A.13) 



2 , 9 1 

- + 21og/i 2 + -- 
e u 



7 9 - 2 = C 2 (N) 

J dxdy 5(1 — x — y)((4 — 6x + 3x 2 ) log(xm,j 2 + (— x + x 2 )g 2 + yA 2 ) + 

(A. 14) 



(2-x) 2 (-x + x 2 )g 2 



xrriq 2 + (—x + x 2 )q 2 + yA" 



•si 



C 2 (G) 



5 5, 2 



J dxdy 5(1 -x - y)((l + 4x - 4x 2 ) log(A 2 + (-x + x 2 )g 2 ) , (A.15) 



2 2, 2 
-3i~3 l0g/i + 



^2 = C 2 (G) 

4 y dxdy 5(1 — x — y) (x(l — x) \og(m g 2 + (— x + x 2 )q 2 ) 
F 9 s = C(N) 



(A.16) 



4 4 , 2 
-3i-3 l0g/i + 



8y dxdy 5(1 — x — y)(x(l — x) log(m g 2 + (— x + x 2 )q 2 ) , (A. 17) 



^4 = C(JV) 



2 2, 2 
"3i-3 l0g/i + 



2 y" dxdy 5(1 -x-y)((l-2x) 2 log(mg 2 + (-x + x 2 )g 2 ) . (A. 18) 

Notice that Fig. ^| (3) does not contribute to the q self energy corrections. 
The expressions for F g \ and F g4: are combinations of three or two diagrams. 



B The Super-Oblique Corrections 

The super-oblique corrections comes from the gluon self energies F g3 , F g4 
(Eq. |A"T7| , Eq. jOg ) and g self energy F g2 (Eq. Then the relative 
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difference of g and g is 

9-9 
9 



a 

Air' 



a F g3 + F t 



ob 



94 



F- n 



<f2 



Arc 



(B.l) 



In the limit m 2 <C q 2 and g 2 <C m^ 2 , these three integral can easily be done 
and give 



F§2 


= C(N) 




= C(N) 




= C(N) 



-- - 2 log /i 2 + 21ogm^ 2 - 1 



4 
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-5= - o^ogfi +-\ogq - — 
Se 6 6 9 



2 

3i 



2 2 
-log/i 2 + -logm/ 



(B.2) 
(B.3) 
(B.4) 



Substituting them into Eq. |B.1| and taking C(N) = 1/2, we have 



#3 ~ #3 
93 



log — 



«3 

12tt 



log- 



(B.5) 



Here we use SU(3) as an example and sum over the contribution from both q^ 
and (Jr. If we want to separate (Jl and (Jr (since their mass can be different), 
then the result is half as big, which is just the the formula given in the paper 
(Eq. |2.3| ) when q 2 ~ m~. The super-oblique corrections we give here take into 
account the finite contribution even from the light quark states. They are 
not included in the renormalization of the coupling constant in an effective 
field theory approach as was done in Ref . |TT| . However, the additional loop 
from light states generates a universal finite term which one might as well 
include in the definition of the oblique parameter, as it will in any case 
be included for the physical predictions (assuming one is incorporating the 
nonlogarithmically enhanced terms). 



C Gluon Emission For q —> qg 

In order to cancel the infrared divergences (which come from the massless 
gluon) in the vertex and self energy corrections, we need to include q — > qgg in 
calculating the decay rate for q — > qg. There are three diagrams contributing 
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to the gluon emission in the process q — > qg (Fig. |7|). The matrix element 
M = Mi + M 2 + M 3 , where 

Mi = -^jf q (p\ ™«hA?(k) J + * 2 + ™^ 2 (l + K)t"v~ 9 (p, m~ g ), (C.l) 
Mi = -^ q (p', m q )(l + 75)t fl (g i 2 ^^ 2 ^(fc)^(p, m g ), (C.2) 

The decay rate for on shell q decay q — > qgg is given by 

dF = (2^32^7 ' M ' 2d,ldS - ( &4 » 

where \M\ 2 can be expressed in terms of s = (p + k) 2 ,t = (p' + k) 2 ,u = 
(p + p') 2 , (s + t + u = nig 2 + nig 2 + m 2 ) rriq, m g , m q and A. Integrating over 
the whole phase space, we can get the correction F emiss i on , which is given by 

r ^ q^qgg ~ rO-^emission; (Cy .5) 

47T 

where To is the tree level decay rate for q — > qg: 

s C 2 (N)g 2 



27im q 2 



IPI 



E q E g + \p\ 2 . (C.6) 



Here |p|, E q , E q are all measured in the rest frame of q. 



D 0(a) Corrections To Couplings In the Ex- 
act SUSY Limit 

In the SUSY limit, when m q = m q and m g = 0, supersymmetry is exactly 
preserved, so the supersymmetric Yukawa coupling ^ is equal to the Stan- 
dard Model gauge coupling gi (i = 1,2,3 indexing U(l), SU(2) and SU(3) 
respectively). Using the formulas given in Appendix A for the vertex and 
self energy diagrams, we show explicitly in Table [| and |p that to the first 

d # in Table | and | means any finite number coming from calculation. 
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C 2 (N) C 2 (G) 





logfi 2 


logA^ 


# 


logfi 2 


logA^ 


# 




1 


2 


5 


-1/2 


-1 


-5/2 




1 





3 


-1/2 





-3/2 


F V 3 











3/2 





5/2 


F V 4 











1/2 





3/2 


F q i 


-1 


-2 


-5 











F,2 


-1 





-3 











l/2Fgi 











5/6 


-5/6 





l/2F g2 











-1/3 


1/3 





-^total 











3/2 


-3/2 






Table 4: 0(a) Correction to the gauge coupling constant g in the SUSY limit 



order, g\ = g\. The vanishing of the C2{N) part is guaranteed by the Ward- 
Identity. The coefficient of log A 2 is the same as log/z 2 , which had to be 
true on dimensional grounds. Notice that there are still infrared divergences 
(the coefficient for log A 2 is not zero) . Presumably these will cancel against 
emission diagrams and other diagrams which don't appear away from the 
super symmetric limit, like the collinear divergence coming from g — > gg. 

E The Cancellation Of The Infrared Diver- 
gences 

Some of the vertices and self energies corrections have infrared divergences 
when the momentum goes to zero. That is why we introduce an infinitesimal 
gluon mass A as a regulator. This infrared divergences cannot be cancelled 
unless we add the contribution of q — ► qgg to q — > qg process, so that the 
decay rate for q — > qg is finite. This is also necessary because experimentally 
we cannot separate these two processes due to the finite resolution of the 
detector. 

First let us study the infrared divergence of the M^M^ [M\ is given in 
Eq. |(J.1| ) in \M\ 2 for q — > qgg. When the gluon momentum k is very small, 
we can find the infrared divergences by keeping only the leading order term 
in k in both the numerator and denominator. Under this approximation, we 
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C 2 {N) 






C 2 (G) 




log/i 2 


log A 2 


# 


log/i 2 


log A 2 # 


F vl 


1 


2 


4 


-1/2 


-1 -2 










2 





-1 


F V 3 











1/2 


-1/2 


F V 4 











2 


-1/2 3 


l/2F ql 


-1/2 


-1 


-5/2 








l/2F q2 


-1/2 





-3/2 








l/2F ql 


-1 





-2 








l/2F q2 


1 


-1 











l/2F g 











-1/2 


1/2 


-Ptotal 











3/2 


-3/2 



Table 5: 0(a) Correction to the Yukawa coupling constant g in the SUSY 
limit 



can treat gluon as massless until we do the final integral over & , which would 
be divergent in the lower limit if the gluon is really massless. We then let 
the gluon have an infinitesinal mass A, which would then provide a lower cut 
off for the integral. For a massless gluon, k = (ko, k) with ko = |k|. We get 

[(p + k) A — m q \ 

= -^wsw ^y . (e.i) 

The 4C 2 (N)g 2 p ■ p' is the matrix element squared for the tree level q — ► 
qg, which will give Tq after we perfom the phase space integral of q and g. 



2 



According to Eq. |C.5| , we have 



-C 2 (N)(4tt) 



a (2vr) 3 2A;o (p 1 ■ k) 2 
o / k^dk r dQ m q 2 
4n 2 k J 4^ (p' ■ k) 2 



C 2 (N)(4n) 2 J x ^ = 2C 2 (N) log A 2 . (E.2) 



'-^emission is tne contribution of the F 0miss i O n from MiM-* 
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Similarly we have 



Emission = 2C 2 (iV)l0gA 2 , (E.3) 

Mission = 2C 2 (G)logA 2 . (E.4) 



The analysis of the cross term M ; Mj* is a little bit more complicated, but 
similar. Considering -fission as an example, we getB 

Mission = (C 2 (N) - l -C 2 {G))W f <lH ' 2 ' r<l 



2 ^ Jjy ' Jx (2 7 r) 3 2fc (p / ■ k)(q- k) 
I^v/^-vnw. n2 [ k 2 dk r dfl m q 2 + m q 2 — m g 2 



(C 2 (N) - -C 2 (G))(Att) 



2 JxA^klJ Art (p> ■ k)(q ■ k) 

1 / n , ,2 f 1 1 1 m Q 2 + m q 2 ~ m g 2 



(C 2 (N) - -C 2 (G))(An) J x _ J q de^T^ 

-2(C 2 (N) - -GAG)) logA 2 f 1 d£ , ( , 

2 Jo t^ z m g z + t^(m q z — m q z — m g z ) + m, 



(E.5) 

where to get from step 2 to step 3, we introduced the Feynman parameter 
£ and performed the dfl integral. In the last step, we keep only the term 
proportional to log A 2 and neglecting all the other terms. Similarly, we get 



F 13 
emission 



-C 2 (G) log A 2 / d£ 
jo 



1 2 2 5 

1 m q — m q — m g 



£ 2 (2m g 2 + 2m q 2 — m q 2 ) + £(m q 2 — 3m q 2 — m g 2 ) + m q 2 ' 

(E.6) 



^emission = —C 2 (G) log A / d£ — 9 — — 9 2 ^ 5" ( E ' 7 ) 

jo £ z m q + t,(m g z — m z — m q z ) + m q z 

Now that we have calculated the infrared divergences in the gluon emis- 
sion, let us study the infrared divergences of the qqg vertex corrections when q 
is on shell. Only F V i, F V3 , -^V4 have infrared divergences. For F V i (Eq. |A.5| ), 
the denominator of the second integral is divergent only when x — > 1 and 
y, z — > 0. In this region, we can set x = 1 and y = z = in the numerator of 

fk = k/k 
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Eq. |A.5| . We can also set x — 1 in the A 2 term in the denominator. Keeping 
only the term proportional to log A 2 and neglecting all the other terms, we 
have[fT7[] 

Fvi = (C 2 {N) - 1/2C 2 (G)) 

dxdy 2( mij 2 + m q 2 -m § 2 ) 

(y — xy)rriq 2 — y(l — x — y)m g 2 + (1 — x)(l — x — y)m 2 + A 2 

(E.8) 

Let y = (1 — x)£, to — (1 — x), this expression becomes 

F V i = -(C 2 (N) - 1/2C 2 (G)) 

d^/W 2 ™fX™\~™f p , 21 ^ x2 - (E.9) 

Jo to 2 [t.mf - 4(1 - 4) m 9 + (1 ~ i) m q \ + A 

Performing the uj 2 integral and let ^ — >• 1 — ^, the result is 

X n \2 f 1 ac m~ 2 + m 2 -m~ 2 



F Y1 = (C 2 (N) - -C 2 (G)) log A 2 / de 
2 Jo 



t^rrig 2 + £ im 2 — rriq 2 — m g 2 ) + rriq 2 ' 

(E.10) 

which is exactly the same as —-^emission/2- 

Performing the same calculation, we see that the infrared divergences 
of 2F V 3 an d 2F V 4 cancel -F 2 mission and ^emission, and the infrared divergences 
of F q i, Fg 2 and F g cancel F^ ission , F 2 mission , F e 3 mission . Thus we have shown 
explicitly that the infrared divergences of the vertices, self energies and gluon 
emission do cancel with each other and give an infrared convergent result at 
the end. 



F The Super-Oblique Corrections From a 
Messenger Sector 

In this section, We consider the super-oblique contribution from a sample 
messenger sector which is in the 5, 5 or 10, 10 representation of SU(5). 
We use mf, m s to stand for the mass of the fermion and the scalar, where 
nif = AS* and m s = \J X 2 S 2 ± AFs. In addition, we define m 2 /m 2 — 1 ± x, 
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thus x = XF S /X 2 S 2 . In the limit q 2 < m 2 , and q 2 < m f 2 , Eq. |7CT0| , Eq. |A~TT 
and Eq. |A.18| give 



FS2 



C(N) 



2 log ji 2 + 2 log m 2 



-m^ + 4m 2 m 2 — 3mf — 2m\ log(m 2 /m 2 



-m 2 + mf) 



2\2 
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C(N) 
C(N) 



4 

3i 

2 

3i 



4 4 
-log/i 2 + -logrn 2 

2 2 
-log/i 2 + - logrn 2 



(F.l) 

(F.2) 
(F.3) 



Substituting them into Eq. [B.l| and multiplying by another 1/2 factor as we 
explained in Appendix O, we get 



9i~9i 



a 
24^ 



1 2 , 3(x 2 + log(l-x 2 )) 
'2 l0g(1 - x) + 4^ 



(F.4) 



The correction is large only when x is very close to 1. In the limit that 
\x\ 1, the correction is 



9i-9i 
9i 



a-. 



1927T 



-X 



;f.5) 



Notice that for U(l) coupling, we need to multiply by the hypercharge 
squared. Also, there is an additional factor of 2 for the U(l) since there is no 
C(N) = 1/2 multiplicity factor. If the messenger sector is in the 5, 5 or 10, 
10 of the SU(5) representation, we can decompose it under SU(3) x SU(2) Y 



5 = (3,l)_ 1/3 + (l,2) 1/2 , 
10 = (3,2) 1/6 + (3,l)_ 2/3 + (l,l)i. 



;f.6) 
;f.7) 



Each SU(3) triplet or SU(2) doublet will contribute to the super-oblique 
corrections. In 5, there is one SU(3) triplet and one SU(2) doublet, while 
for 10, there are three SU(3) triplets and three SU(2) doublets. Thus the 
multiplicity factor is 2N, where N = (n 5 + 3nio) if we have n 5 pairs of 5 and 
5, nio pairs of 10 and 10. The super-oblique correction from the messenger 
sector is then 



9i-9i 
9i 



Oil 

96^ J 



Nx 2 



(F.8) 
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Substituting the actual numbers into Eq. [F.8|, we get 



SU(3) - — — = 0M%Nx 2 , (F.9) 

93 

SU(2) 32^2 = QM1%Nx 2 } (F.10) 

U(l) g l^ll = 0.0056%iVx 2 , (F.ll) 
9i 

which is very small. 
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Figure Captions 



Figure 1. One-loop diagrams for the gqq vertex. The arrow shows the glow 
of quark number, q can be either q L or q R . F V ; corresponds to the 
correction from diagram (i), i = 1, 2, 3, 4. Similar notations are used in 



the diagrams below. 




Figure 


2. 


One-loop diagrams 


for the qqg vertex. 


Figure 


3. 


One- loop diagrams 


for g self energies. 


Figure 


4. 


One-loop diagrams 


for q self energies. 


Figure 


5. 


One- loop diagrams 


for q self energies. 


Figure 


6. 


One- loop diagrams 


for g self energies. 


Figure 


7. 


Feynman diagrams 


for q^ qgg. 


Figure 


8. 


Number of events per year for e + e - 



and luminosity L = 5 x 10 33 cm 2 s 1 



tig at E cm = 1000 GeV 



Figure 9. Number of events per year for e + e — > bbg at E cm = 1000 GeV 
and luminosity L = 5 x 10 33 cm _2 s _1 . 

Figure 10. SU(3) non-oblique corrections for on shell squark decay q — > qg. 
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(1) 



(2) 



(3) 



(4) 



Figure 1: One-loop diagrams for the gqq vertex. The arrow shows the glow 
of quark number, q can be either q\ or q&. Fy; corresponds to the correction 
from diagram (i), i = 1,2,3,4. Similar notations are used in the diagrams 



below. 




(1) (2) (3) 



Figure 2: One- loop diagrams for the qqg vertex. 



q 




Figure 3: One-loop diagrams for g self energies. 




(1) (2) 



Figure 4: One-loop diagrams for q self energies. 
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Figure 10: SU(3) non-oblique corrections for on shell squark decay q — > qg. 
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